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Abstract. - We study the scaling regimes for the Kardar-Parisi-Zhang equation with noise 
correlator R(q) oc (1 + wq~ 2p ) in Fourier space, as a function of p and the spatial dimension 
d. By means of a stochastic Cole-Hopf transformation, the critical and correction-to-scaling 
exponents at the roughening transition are determined to all orders in a (d — d c ) expansion. We 
also argue that there is a intriguing possibility that the rough phases above and below the lower 
critical dimension d c — 2(1 + p) are genuinely different which could lead to a re-interpretation 
of results in the literature. 



The Kardar-Parisi-Zhang (KPZ) equation, which was originally introduced to describe 
growth of rough surfaces [Q, displays generic scale invariance, as well as a non-equilibrium 
roughening transition separating a smooth from a rough phase above the lower critical dimen- 
sion. As a consequence of its mapping to the noisy Burgers equation to the statistical 
mechanics of a directed polymer in a random environment ||, as well as to other interesting 
equilibrium and non-equilibrium systems (for recent reviews, see Refs. HH), the KPZ problem 
has emerged as one the fundamental theoretical models defining possible universality classes 
for non-equilibrium scaling phenomena and phase transitions. 

In one dimension, the roughness and dynamic exponents, x an( i z -> have long been deter- 
mined exactly by means of the dynamic renormalization group (RG), utilizing the symmetries 
of the problem Q . Furthermore, it has been demonstrated that the associated scaling functions 
can be computed to high precision by means of the self-consistent one-loop or mode-coupling 
approximation || fjj. For d > d c , a two- loop RG calculation || indicated that the critical 
behavior of the roughening transition might be described by an exact set of exponents as 
suggested earlier on the basis of scaling arguments B. Using a directed-polymer repre- 
sentation, Lassig was able to demonstrate the validity of this statement to all orders in a 
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(2 + e) expansion |Io| (see also Refs. jy], [l2|]). In addition, these results were evidence for a 
upper critical dimension d uc = 4 of the roughening transition, and suggested that the ensuing 
strong-coupling rough phase was not accessible within perturbation theory. This is in accord 
with the divergence of the coupling constant g upon approaching the lower critical dimension 
from below H|. 

The scaling behavior in the strong-coupling rough phase above d c has been very controver- 
sial. Based on very different assumptions and analytic approaches, diverse values for the scaling 
exponents have been postulated, see, e.g., Refs. II £§. In addition, some authors have 
claimed d uc = 4 to be the upper critical dimension for the rough phase as well [ 16 , |7l |1| , where 
the scaling exponents assume the values known in infinite dimensions |l9|, 20 1. In contrast, 
numerical studies observed merely continuously varying exponents as d was increased j2lj . In 
addition, the validity of the continuum Langevin description has been questioned in this regime, 
and a conceivable breakdown of universality has been conjectured (see, e.g., Refs. p^, |23[]). 

In order to shed light on some of these open issues, we shall find it useful to add a long-range 
power-law contribution to the usual spatially local stochastic noise of the KPZ equation, as 
first introduced by Medina et al. p|. More generally, we investigate the Langevin equation, 

a tS (x, t)=D V 2 s(x, t) + ^ [Vs(x, t)] 2 + C(x, t) , (1) 

with Gaussian noise characterized by zero mean and variance (C(x, t) C( x '> 0) = 2i?o(x — 
x') S(t — t'), which we assume to be local in time, but which may contain spatially long-range 
power law contributions of the form i?o( x — x ') oc |x — x'| 2p ~ d . Thus, typically we have 
-Ro(q) = D (l + w |q| -2 '') in Fourier space. Notice that setting u = — Vs in Eq. (||) leads to 
the noisy Burgers equation, with the usual locally conserved noise for w = or p = 0, but with 
non-conserved noise for p = 1 (termed model B in Ref. Q]). We find as the most important 
effects of adding the long-range noise term: (i) the lower critical dimension for the roughening 
transition is shifted upwards to d c = 2(l+p); (ii) above d c , there are two subtly distinct 
regimes for the smooth phase, characterized by different correction-to-scaling exponents; (iii) 
below d c , there are two distinct rough regimes, governed by the short-range and long-range 
noise RG fixed points, respectively, and separated by a line p*(d) in the (p, d) plane. 

In the following, above d c , we derive an exact integral equation for the noise correlation 
function by means of a stochastic Cole-Hopf transformation. Based on the ensuing minimally 
renormalized RG flow functions, we compute critical exponents at the roughening transition to 
all orders in a (d—d c ) expansion, and determine the exact scaling exponents in the smooth phase 
above d c . We demonstrate that the strong-coupling rough phase above d c is perturbationally 
inaccessible, but probably characterized by a spatially local noise correlator. Below the lower 
critical dimension d Cl we determine the scaling exponents at the long-range noise fixed point 
exactly, provided such a non-trivial finite fixed point exists. We discuss different analytical 
RG approaches to such a strong-coupling fixed point, and obtain an approximate expression 
for the separatrix p*(d). Finally, we discuss whether the rough phases above and below d c , are 
"continuously" connected as a function of space dimension d and the magnitude of the noise 
correlation exponent p. We shall argue that there is a intriguing possibility that these two 
rough phases are genuinely different. 

A convenient starting point for a systematic analysis of a Langevin equation like Eq. ([l]) is 
its reformulation in terms of a dynamic generating functional § , 



J[s,s] = d d x dt s(x,t) 



d t s(x,t)-DS7 2 s(x,t) (2) 



[Vs(x, tf - J d d x'R ( X - x') ,5(x', t) 
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Fig. 1. - Bethe-Salpeter equation for the noise correlator. 



which allows the calculation of expectation values by means of path integrals with the statis- 
tical weight exp(— J[s, s]). Upon directly applying a Cole-Hopf transformation to the KPZ 
equation, Eq. (Q), one obtains a diffusion equation subject to multiplicative noise which is 
often interpreted as a directed polymer in a random potential || . If we recast the same idea in 
terms of the above dynamic functional the corresponding stochastic Cole-Hopf transformation 
[ fi"T| — in the appropriately discretized version of Eq. (^) in the Ito representation — reads 



n(x,i) = - exp{| [s(x,i) +L>i? (0)t]} , 
n(x,i + T) = s(x,t + r)exp j-|[s(x,t) +£>i? (O)t]j . 



(3) 



This leads to a dynamic generating functional, 



J[h,n}= d d x dt\h(x,t) 9 t n(x, t) - D V 2 n(x, t) 



Dg 2 



(4) 



d d x' n(x, t) n(x, t) R (x - x') n(x', t) n(x', t) 



very reminiscent of the field theory for diffusion- limited pair annihilation [|25| . A remarkable 
feature of Jj[n, n] is that there exist no loop diagrams contributing to a renormalization of the 
diffusion propagator. Hence the dynamic exponent is z — 2 whenever standard perturbation 
theory is applicable. This leaves us with the renormalization of the four-point noise vertex, 
which is formally achieved to all orders in the perturbation expansion via a Bethe-Salpeter 
equation, as graphically depicted in Fig. |l]. Analytically, this leads to an integral equation for 
the renormalized noise correlator R, 



i?(k,k' ;A i 2 ) = i? (k-k') 



4 



d d P flo(k-p) 



(27T) d p 2 + H 2 



i?( P ,kV) 



(5) 



where /i 2 = iu/2D + q 2 /4, and q and lo are the momentum and frequency transfers from right 
to left in the vertices of Fig. [I], while k — k' denotes the momentum transfer from bottom 
to top. The standard perturbation expansion in terms of g 2 is equivalent to the Neumann 
series for this Fredholm integral equation. The exact relation rtq) can be used to determine the 
scaling function for the renormalized noise self-consistently p6| . 

In this letter, we focus on the asymptotic scaling exponents. Upon inserting the noise 
correlator into the Bethe-Salpeter equation (^), one readily obtains the full renormalization 
of the short-range noise strength D within the minimal-subtraction procedure, which here 
leads to a formal double expansion with respect to e = d — 2 and p. On the other hand, the 
non-analytic long-range contribution itself is not renormalized perturbationally. In terms of 
the renormalized counterparts u and v for the bare couplings g 2 and wg 2 , respectively, the 
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ensuing exact RG (3 functions are actually precisely those of the one-loop theory, 

I3 U = eu- [u + vf/2 , v = (e-2p)v. (6) 

For d > d c = 2(1 + p), Eq. (||) shows that v — > asymptotically. Therefore, the usual 
short-range noise KPZ scenario applies [I], [24|, ^]. 

In addition to the trivial fixed point u* = and the strong-coupling fixed point it* = oo, 
there appears an unstable critical fixed point (u c ,v c ) = (2e,0), marking the location of a 
non-equilibrium roughening transition above the lower critical dimension d c . At this O(s) 
fixed point, the general scaling relation x + z = 2 is valid, and because of z c = 2 we find 
the marginal roughness exponent Xc = |j, 0. Furthermore, the negative eigenvalue of 
the associated stability matrix defines the crossover or inverse correlation length exponent 
4> c = v^ 1 = d — 2, whereas its positive eigenvalue yields the correction-to-scaling exponent 
uj c = d — d c . Assuming that the Chayes-Fisher bound v > 2/d |^7j applies for the crossover 
length scale in this problem (equivalent to a directed polymer in a random environment), one 
finds that this lower bound is reached in d uc — 4 dimensions, which therefore constitutes the 
upper critical dimension for the roughening transition, beyond which <p c = 2, v c = 1/2 fi"o| , psf . 

In the smooth phase, both couplings mO and v — > 0, and a more careful analysis of the 
RG flow of the ratio w — v/u is required (details will be presented elsewhere [^6)). One finds 
that actually w — > oo, and thus the smooth phase, as opposed to the roughening transition, is 
characterized by the long-range algebraic noise, with roughness exponent Xsm = l+p—d/2 < 0, 
and z sm — 2. Furthermore, two distinct scaling regimes emerge, distinguished by different 
correction-to-scaling exponents: For d c < d < 2(1 + 2p), one finds io\ = d — 2(1 + p) and 
uj2 = 2(1 + 2p) — d, whereas for d > 2(1 + 2p), lo\ = d — 2(1 + 2p) and LO2 = 2p. Precisely at 
d = d c , i.e., for p — (d — 2)/2, there appears a fixed line, which is unstable for w < 1, and 
stable for w > 1. 

Finally, the RG analysis also tells us that the rough phase emerging above u c is a genuine 
strong- coupling phase in the sense that it remains inaccessible to perturbative methods, even to 
all orders in e. Numerical solutions of the flow equations also suggest that w — > in the rough 
phase; hence the algebraic noise correlations appear to be irrelevant in the strong-coupling 
regime. 

Up to now we have restricted our discussion to dimensions larger than the lower critical 
dimension, where there is a kinetic roughening transition from a smooth to a rough phase. 
If one tries to extend the above analysis to d < d c — 2(1 + p), where there exists no 
roughening transition, Eq. (0) implies that v — > 00 at long length scales. Thus, the minimally 
renormalized perturbation theory based on the stochastic Cole-Hopf transformation breaks 
down, and one must resort to other methods. Fortunately, however, we can draw some 
important exact conclusions already from the general structure of the field theory As 
a consequence of (a) the underlying Galilean invariance, which fixes the renormalization of 
g, (Jo) the fact that the non-analytic noise term proportional to Dw cannot renormalize, and 
(c) the momentum dependence of the non-linear vertex, there are merely two independent 
renormalization constants to be determined, namely for the renormalized fields, sr = Z 1 ^ 2 s 
and the renormalized diffusion constant, Dr = ZpD. Upon defining 7 = /i<9^|oln Z and 
C = /i9 M |oln^D, the RG (3 functions can be expressed as f3 u = (d — 2 — 7 — 2Q)u and 
Pv = (d — 2 — 2p — 3C)w. Then by solving the RG equation for the correlation function 
near a stable RG fixed point, we furthermore identify x = (2 — d + 7*)/2 and z = 2 + £*. The 
existence of a non-zero, finite RG fixed point u* then immediately leads to the scaling relation 
X + z = 2 Similarly, at any long-range fixed point < < 00, is fixed, giving the exact 
values 



zi r = (4 + d-2p)/3 , X ix = (2-d+2p)/3 



(7) 
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Fig. 2. - Scaling regimes (left) and sketches of the RG flow of the coupling constant g (right) above 
and below d c (p). The dashed-dotted line indicates a approximate separatrix between the domain of 
attraction of the short-range and long-range noise fixed point. 



for the scaling exponents at the long-range fixed point, provided < u* < oo as well. 

In general, these two scaling fixed points compete, and the short-range fixed point must 
evidently be dominant, if z aT < z\ r (and vice versa), which indeed implies (3 V > and hence 
v — > 0. In one dimension, and for purely local noise (w = 0), the non- linear reversible force 
term proportional to g can be shown to fulfill the detailed-balance conditions, guaranteeing 
that the stationary probability distribution becomes V s t oc exp (— J dx[Vs(x)] 2 /2) , as for the 
linear equation. In this Gaussian static theory, there can be no field renormalization (Z = 1), 
hence £ = 0, leading to the familiar one-dimensional KPZ short-range scaling exponents 
Xsr — 1/2 and z sx = 3/2 [Q, [ij. Comparing the latter with the long-range dynamic exponent 
zi r = (5 — 2p)/3, we find that the short-range fixed point remains stable, provided p < 1/4 
[ p4| . After some controversy in the literature, this result has been confirmed by simulations 
for the noisy Burgers equation (29). 

Notice that a minimal renormalization scheme can never arrive at the exact one-dimensional 
exponents x sr = 1/2 and z sr = 3/2. Instead, in order to address the rough phase below 
d c , one may devise a non-minimal renormalization procedure at fixed dimension d and p 
[|| ^6). Alternatively, one may utilize the mapping to the Burgers equation and hence to 
driven diffusive systems, for which a well-defined (2 — d) expansion exists. Adding long-range 
correlated noise, this actually leads to the identical stability condition p < 1/4 for the short- 
range fixed point in d = 1 ]3Q] . In the long-range regime, the case p = 1, corresponding to 
the Burgers equation with non- conserved noise, is accessible through an e expansion below the 
upper critical dimension d uc = 4 of this model The dynamic exponent here is actually 
obtained to all orders in e, and reads z\ Y = (2 + d)/3. 

In order to further discuss the implications of the above exact results for the KPZ equation 
with long-range correlated noise it is instructive to consider the scaling regimes in the (d, p) 
landscape; see Fig. ^|. One should notice that there are two qualitatively distinct regions 
separated by the lower critical dimension d c (p) = 2(1 + p). For d < d c {p) there is only a rough 
phase and the short-range (SR) and long-range (LR) noise fixed points compete. Above the 
lower critical dimension there appears a phase transition from a smooth to a rough phase and 
the RG analysis indicates that LR noise constitutes an irrelevant perturbation in the rough 
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phase. 

From the above analysis we know the values of the scaling exponents in the domain of 
attraction of the LR fixed point exactly. Since the scaling exponents at the SR fixed point are 
independent of p, the only additional information needed to determine the SR exponents is the 
separatrix p*(d) between the domains of attraction of the LR and SR fixed points. Of course, 
for the latter to be true one has to require that the exponents are continuous at p*(d). This 
has been shown explicitly for d = 1, and seems very reasonable in general. From the analysis 
in this letter we cannot determine the exact form of the separatrix. Yet, we actually can locate 
some points on this curve. In d = 1, the short-range fixed point is stable for p < p* = 1/4; 
for p=l (the Burgers equation with non-conserved noise), and below d uc = 4, there is 
only the long-range regime. If we assume that the separatrix p*(d) between the short-range 
and long-range regimes extends up to four dimensions, then it definitely contains the points 
(1, 4), (1, 1/4), and probably also (0, 0) in the (p, d) plane. A simple linear interpolation yields 
the function p*(d) ~ d/A (dashed line in Fig. |J); recent computer simulations in fact found 
p*(2) « 1/2 ]3lj |. Inserting the linear interpolation approximation (which we do not expect to 
be exact), p*(d) « d/A, into Eq. (Jfy yields 

z sr «(8 + d)/6, « (4 - d)/6 . (8) 

These values remarkably coincide with Halpin-Healy's results obtained in a functional RG (l3) 
and also with a more recent perturbative mode-coupling study by Bhattacharjee ||T8"|| . 

There are now two possible scenarios. The first one is that all the scaling exponents for 
the strong-coupling phase are continuous over the whole (p, d) plane and in particular upon 
crossing the line marking the lower critical dimension. This would necessarily imply (if not 
prove) that d uc = 4 is the upper critical dimension for the rough phase in agreement with some 
recent speculations Jl7| but disagreement with computer simulations [ pl| . However, one should 
probably expect to encounter singularities as the lower critical dimension d c (p) is crossed. A 
second more intriguing scenario is therefore that there are two fundamentally different rough 
phases, one below and one above the lower critical dimension d c {p) = 2(1 + p), which we 
call type-I (SR-I) and type-II (SR-II). This would then allow for the scaling exponents to be 
discontinuous at d c (p). Below d c {p) one would have scaling exponents close to those obtained 
from the linear interpolation of the separatrix in Eq. (^). In particular the upper critical 
dimension of the type-I rough phase is d uc — 4. Above the lower critical dimension d c {p) one 
would have a different set of exponents (e.g. those given by the numerical simulations) and 
the upper critical dimension must not necessarily be equal to 4 or any other finite value. At 
present there is proof for neither of the above scenarios but the following observations are 
quite indicative. First, from the RG analysis exploiting the Cole-Hopf transformation we have 
learned that the strong-coupling phase above d c (p) is not accessible by perturbation theory 
even to infinite order. On the other hand the rough phase at d — 1 is accessible by standard 
perturbation theory using a mapping of the KPZ equation to a driven diffusion model [ j30f . 
Second, for p = an explicit two-loop calculation || shows that the fixed point value of the 
coupling constant g approaches infinity as the lower critical dimension approaches 2 from below. 
Third, this scenario would provide a coherent picture for most of the available numerical and 
analytical results for the KPZ equation. Some of the analytical approaches (e.g. functional 
RG and mode-coupling theory) are self-consistent theories and hence necessarily start out 
with correlated noise. This in effect shifts the lower critical dimension upwards, and it is 
well possible that this automatically constrains the results to the SR-I phase as opposed to the 
SR-II phase these studies were supposedly aiming at. It is difficult to judge on other analytical 
methods and how they fit into the scheme discussed here. Lastly, there have been suggestions 
for a breakdown of a continuum description, and perhaps even universality, in the rough phase 
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[ p2| |2^| . Physically, in a microscopically rough regime, the underlying lattice as well as details 
of the dynamic rules may well be important; in addition, one might question the existence of 
a well-defined short-distance expansion |HJ in this phase. Clearly, these issues require further 
clarification through approaches that extend beyond the continuum equation, Eq. ([!]), and 
equivalent field theory methods. 
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